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For our English DUG-Members : 

One of the results of the DUG-Meeting in Nottingham was the wish of some 
members to have more contacts within England. Two participants of this mee- 
ting are willing to build up a regional DERIVE Group to exchange informati- 
ons and experiences within the UK. I support these efforts and hope to ex- 
tend the relations to my British friends. I want to thank all participants 
of the meeting in Nottingham for their co-operation and especially for 
their cordiality. These days in England were a very special event for my 
daughter , my wife and me. We will have best memories on this meeting. 


DERIVE Group South England, please contact 

Kevin Williamson 

South East Sixth Form College 

Runnymede Chase 

Benfllet, Essex SS7 1TW 

DERIVE Group North England, please contact 

Terence A. Etchells 
98 High Street 
Belmont, Bolton BL7 8AL 

Best greetings to Kevin and Terence. I wish you both much success and I’m 
sure to receive many interesting contributions for the DNL in the future. 

Maybe that there are some DUG-Members in other countries who want to found 
a regional DERIVE-Group . We would be glad to encourage and to support your 
ideas . 


A short report of the International Spring School 

The first International Spring School on the Didactics of Computer Alge- 
bra took place in Krems, Austria, from April 27 - 30, 1992. It was an in- 
spiring and innovative meeting of 40 experts. But there was a tragic event, 
too. Dr Dennis Tanner, London, had a bad accident on the first evening from 
which he couldn't recover. So I'm sorry to tell you, that our friend Dennis 
died 10 days later in a hospital. I've met Dennis twice for some hours and 
I will never forget his humour, his spirit and his vitality. 

In the next DNL you will find more details about the Krems Meeting. The 
proceeding is to be published (about 320 pages! !) . On the last page you can 
find a list of the participants, who gave a lecture. If anybody wants to 
receive some details about one or another item, then please inform the edi- 
tor . 


I must thank all participants for their friendship. I've got a lot of 
letters since then and I think that this week was a very special week in 
our life. I can't describe it really, but there was more than a common in- 
terest in DERIVE. We all felt, that something more happened in Krems. Best 
greetings to all of you and I am sure to meet you again somewhere. 

My special thank to Dr. Heugl and Dr. Kutzler. Your organization was 
perfect. I'll enjoy the next Spring- or Summer School. 
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Lieber Derive Anwender! 


Dear Derive User, 


Ich mochte mich dafiir entschuldigen, dass dieser 
DNL so spat im Juni erscheint. Aber es gab so viel 
zu tun in letzter Zeit. Vor allem die Aufarbeitung 
der DERIVE Spring School in Krems und die Zu- 
sammenstellung des Tagungsberichtes hat einige 
Arbeit gemacht. Wenn ich Ihnen sagen kann, in 
welcher Form dieser Bericht erscheint, werde ich 
Sie natiirlich hier informieren. 

Der angekiindigte Artikel iiber italienische 
DERIVE-Aktivitaten musste fur den nachsten DNL 
aufgeschoben werden. Erstens war die Ubersetzung 
fur mich etwas milhsam: italienisch - deutsch - 
englisch und zweitens liegt mit H. Appels Wunschli- 
ste bereits ein didaktischer Aufsatz vor. Und ich bin 
der Meinung, dass der DNL keine nur Mathematik- 
didaktische Schrift werden soil, obwohl didaktische 
Themen nach wie vor willkommen sind. Daher 
auch an dieser Stelle meine wiederholte Aufforde- 
rung an alle DERIVE-Anwender: bitte berichten 
Sie von Ihrer ganz speziellen DERIVE-Anwendung ! 

Unmittelbar nach Fertigstellung des DNL#6 
erhielt ich aus Hawai eine Diskette mit den Anfra- 
gen und Antworten aus dem DERIVE Electronic 
Bulletin Board von September 91 bis Mai 92. Zu- 
gleich erhielt ich die Zusage, monatlich mit den 
neuesten Informationen versorgt zu werden. Vielen 
Dank, mahalo nui loa! 

Juni ist Rosenmonat, daher habe ich einige Mau- 
rer-Rosen in diesen DNL eingestreut (P.Maurer, 
"A rose is a rose ... " Amer. Math. Monthly, 94). Im 
nachsten DNL konnen Sie die DERIVE Wurzeln 
und Zweige - und Dornen— dazu finden. 

Ich wiinsche alien DERIVE-Freunden einen scho- 
nen Sommer und freue mich bereits jetzt auf ein 
„ Wiederlesen “ im Herbst. 


I have to apologize for publishing this DNL so 
late in June. But there was so much to do last time, 
above all the working up of the DERIVE Spring 
School held in Krems and then collecting and edi- 
ting the proceedings. I will let you know in the DNL 
when and in which form this book will be available. 

The announced contribution concerning Italian 
DERIVE activities from R.M.Castelletti had to be 
postponed for the next issue because of two rea- 
sons: at first it was a little bit difficult for me to 
translate Italian — » German —> English and on the 
other hand there is already a didactic article in this 
DNL. In my opinion the DNL should not become an 
only maths didactic magazine, although didactic 
items are welcomed now as before. Therefore once 
more my invitation to all DERIVE users: please 
report about your special DERIVE-application(s)! 

...When DNL#6 had just been finished I received a 
diskette from Hawaii containing requests, answers 
and comments from the DERIVE Electronic Bulle- 
tin Board from Sepember 91 through Mai 92. I got 
the promise to be informed about the latest news 
from the BBS monthly from now on. Many thanks 
Hawaii, mahalo nui loa! 

June is the month of roses, therefor I've strewed 
some Maurer roses for you in this DNL (P.Maurer, 
"A rose is a rose...”, Amer. Math. Monthly, 94). In 
the next DNL you will find the roots and branches - 
and thorns. 

I wish a fine summer to you and your families 
and I'm looking forward to a "reading again" in 
autumn 

With my best regards 



Mit den besten Griifien 





EDITORIAL 
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The Derive-News-Letter is the Bulletin of 
the Derive User Group. It is published at 
least three times a year with a contents of 
30 pages minimum. The goals of the D-N- 
L are to enable the exchange of experi- 
ences made with Derive as well as to cre- 
ate a group to discuss the possibilities of 
new methodical and didactic manners in 
teaching Mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 


Contributions 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE - Newsletter will be. 


Preview: Contributions waiting to be published 

Logic with DERIVE 
Maurer Roses 
Probability Theory 
DERIVE in Italy 


will be published in September 1992 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, DXust 1, AUSTRIA 
Richtung: Informationsblatt der DERIVE User Group 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Prof. Dr. V. WeiRpfennig, Passau 

During a course on polynomial factorization that I gave at the University of Passau this winter, I tried 
to use DERIVE for computing some examples. In this connection, DERIVE turned out to be of very 
limited use only, since most factorization algorithms work by modular techniques. On the other hand, 
DERIVE does not admit polynomial calculations (such as GCD calculations or various factorizations) 
over integers modulo a prime p. 


The biggest surprise, however, was the polynomial 

J{x) = (x 3 + 3x 2 - l)(x 3 - x + 1) = x 6 + 3x 5 - x 4 - 3x 3 + 3x 2 + x - 1 

that I chose as a convenient example for factorization in Z[x] using Berlekamp's method over Z/2 and 
Hensel lifting. When I asked DERIVE to factor /(*) over the rationals (or even over the complex num- 
bers), DERIVE insisted that/(*) is irreducible! 


Another example of this kind is 

g(x)=(x 4 - x 2 + 5x - l)(x 3 + 4x + 1) = x 7 + 3x 5 + 6x 4 - 5x 3 + 19x 2 + x - 1 
Has anyone else noticed this bug? Probably, these examples are just the tip of an iceberg. 


D-N-L: David Stoutemyer confirmed your question. The algorithms mentioned by you are not imple- 
mented in DERIVE. Maybe future versions will be able to do better factorization. 


As you can see , this problem has been resolved in the meanwhile, Josef 2005: 
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Volker Neurath, Velbert 

... Ich betreibe unter dem Betriebssystem DR-DOS 5.0 die Version 2.07 von DERIVE. 

Das Problem ist, dass bei mir der Befehl Transfer Print Printer nicht funktioniert. DR-DOS bricht 
den Vorgang mit einer Fehlermeldung ab. Zur Zeit drucke ich daher mit Transfer Print File in ein 
File, dessen Inhalt ich fiber Options Execute mit Hilfe einer Batchdatei drucken lasse. 

D-N-L: The reason is an incompatibility of DR-DOS 5.0. Please ask DR for a patch. You will get it 
and then DERIVE will print its files. 

Some weeks later I received another mail from Mr. Neurath: 

...Danke auch fur Ihre Information beziiglich DERIVE/DR-DOS 5.0. Allerdings hat sich dieses Pro- 
blem inzwischen erledigt. Teilen Sie alien Anwendem mit: DERIVE und DR-DOS 6.0 verstehen sich 
in jeder Beziehung prima. Ich setze jetzt seit circa 4 Wochen DR-DOS 6.0 ein und bin bisher sehr 
zufrieden. 
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Arney, Exploring Calculus with DERIVE , Addison-Wesley, ISBN 0-201-52839-8 

DeMarois, College Algebra Laboratories Using Derive, MathWare, 604 East Mumford Drive, 

Urbana, IL, 61801, USA 

Harper, Woof & Hodgkinson, A Guide to Computer Algebra Systems , John Wiley, 

ISBN 0-471-92910-7 

Johnson & Evans, Discovering Calculus with Derive , John Wiley, ISBN 0-471-55155-4 

Salter & Gilligan, Linear Algebra Experiments Using the Derive Program , 

Gilmar Publishing Company, P.O. Box 6376, Cincinatti, OH, 45206, USA 

Wasen & Sjostrand, Svensk Handledning med Ovningsuppgifter och Introduktion till Algebraisk 

Modellering och Datoralgebra: Derive , GraftiTexData, Storgatan 11, 590 40 Kisa, Sweden 

Marlewski, Derive Pomocnik matematyczny wersja 2.0 , Wydawnictwo NAKOM, Poznan, Poland. 

ISBN 83-85060-49-9;A polish manual, with a lot of examples (355 p.) 

Scheu, Arbeitsbuch Computer-Algebra mit DERIVE, Dummler, Bonn, ISBN 3-427-45721-4 
Beispiele, Algorithmen und Aufgaben aus der Schulmathematik 
aus dem Inhalt: Kurvenscharen, Newtonverfahren, Regula falsi, 

Grenzwerte, Programmieren mit DERIVE, Experimente mit 
Primzahlen, .... 

Scheu, Entdeckungen in der Menge der Primzahlen mit DERIVE, PM 34, (Heft 3) 1992 

Hischer (Hrsg.), Mathematikunterricht im Umbruch? Erorterungen zur moglichen 

„ Trivialisierung “ von mathematischen Gebieten durch Hardware und Software ; 

Proceedings, Verlag Franzbecker, Hildesheim 1992. 

A lot of lectures concerning DERIVE (in German). 
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This time we open Mr Setif's Treasure Box and there we can find some jewels glittering like Fractals. 

APV(u , v) was Mr Setif's function to append vectors. Later APPEND(u , v) was implemented in 
DERIVE. So one can replace all APV-appearances by APPEND. 

APV(u , v) :=VECTOR(IF(m_<DIMENSION(u) , ELEMENT(u ,m_) , 

ELEMENT(v,m_-DIMENSION(u))) , m_, DIMENSION (u)+DIMENSION(v)) 

FL(x,y,l ,a,m,n) :=IF(n<0, [] , 

IF(n, [[x+1 -COS(m-7r/3) ,y+l • SIN(m-7r/3) ] ] , 
APV(FL(x,y,l/3,m,n-l) , 

APV(FL(x+l •COS(m-7T/3)/3,y+l •SIN(m-7r/3)/3,l/3,m+l,n-l) , 
APV(FL(x+l •(3-COS(m-7r/3)-V3-SIN(m-7r/3))/6, 

y+1 • (3-SIN(m-7T/3)+73-COS(m-7r/3))/6,l/3,m-l,n-l) , 
FL(x+2-l -COS(m-7T/3)/3,y+2-l -SIN(m-7r/3)/3,l/3,m,n-l)))))) 

FLK(x , y , 1 , m,n) :=APV([[x,y]] ,FL(x,y,l ,m,n)) 

FLA(x,y,l , a,m,n) := IF(n<0,[], 

IF(n, [[x+1 • COS (a+m -77/3) ,y+l -SIN(a+m-7r/3)]] , 
APV(FLA(x,y,l/3,a,m,n-l) , 

APV(FLA(x+l -COS(a+m-7r/3)/3,y+l -SIN(a+m-7T/3)/3,l/3,a,m+l,n-l) , 
APV(FLA(x+l • (3-COS(a+m-7T/3)-V3-SIN(a+m-7T/3))/6, 

y+1 • (3-SIN(a+m-7T/3)+V3-COS(a+m-7T/3))/6,l/3,a,m-l,n-l) , 
LA(x+2 • 1 • COS (a+m • 7 t/ 3) /3 , y+2 • 1 • SIN (a+m - 7 t/ 3) /3 ,1 /3 ,a,m, n-1) ) ) ) ) ) 

FLKA(x,y,l , a,m , n) :=APV([ [x , y] ] , FLA(x , y , 1 ,a,m,n)) 

FLAKE_R(x,y,l , n) :=[FLKA(x ,y , 1 , 0 , 0 , n) , FLKA(x+l ,y,l , 2 -7 t/ 3 , 0 , n) , FLKA(x+l /2 , 
y+1 -73/2,1 , 4 - 7r/3 , 0 , n) ] 


FLAKE_S(x,y,l ,n) := [ FLKA(x , y , 1 ,7r/3 , 0 , n) , FLKA(x+l /2 , y+1 -73/2,1 ,-7r/3,0,n) , 

FLKA(x+l ,y,l ,7T,0,n)] 


Picture 1 (for the DOS-version) 


FLK(-3,0,6,0,1) 

FLK(-3,0,6,0,2) 

FLK(-3,0,6,0,3) 

FLK(-3,0,6,0,4) 


Picture 1 (for DERIVE 6) 


FLK(-6 . 5 ,1 ,6 , 0 , 1) 
FLK(0 . 5 , 1, 6, 0 , 2) 
FLK(-6 . 5 , -3 , 6 , 0 , 3) 
FLK(0 . 5 ,-3 ,6 , 0 ,4) 



Picture 2 
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Picture 3 
FLAKE_R(-1 , 

-1, 

3, 

1) 

FLAKE_R(-1 , 

-1, 
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2) 

FLAKE_R(-1 , 
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First we calculated some iterations of the Koch-curve. The next file shows fractal candelabers: 


CAND(x ,y ,1 , n) := IF(n,[x, y + 2-1; x, y] , 

APV([x , y + 1 ; x+1, y + 1], 

CAND(x +1, y + 1, 1/2, n - 1) , [ [x — 1 , y + 1]], 
CAND(x - 1, y + 1, 1/2, n - 1), [x, y + 1; x, y])) 


ml = multiply length branch 
me = multiply length candle 

CANDMM(x,y,l ,ml ,mc,n) := IF(n,[x, y + mc-1; x, y], 

APPEND([x, y + 1; x+1, y+1], 

CANDMM(x +1, y+1, ml-1, ml, me, n - 1) , 
[[x - 1 , y + 1 ]] , 

CANDMM(x - 1, y+1, ml-1, ml, me, n - 1) , 
[x, y + 1 ; x, y])) 

picture 4 

[CAND(-6, -1, 1, 1), CAND(-2, -1, 1, 2), CAND(2, -1, 1, 3), CAND(6, -1, 1, 4)] 


LMJ [MJ 4 4^ 4^ ^ 
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Picture 5 


[CANDMMC-9, -2, 1, 0.75, 3, 1), CANDMM(-3 , -2, 1, 0.75, 3, 2), 
CANDMM(3 , -2, 1, 0.75, 3, 3), CANDMM(9, -2, 1, 0.75, 3, 4)] 



Picture 6 

CANDMMCO, -2, 1, 0.25, 3, 3) 
CANDMMCO, -2, 1, 0.25, 10, 3) 
CANDMM(0, -2, 1, 0.4, 10, 3) 


hi. 


-b - 1 


hi. 


-b - 1 




CANDMMCO, -2, 1, 0.47, 15, 6) 
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Trickfilme mit DERIVE 
Graphisches Differenzieren 

Hans-Jiirgen Kayser, Diisseldorf 

DERIVE ist hervorragend zur Erstellung veranschaulichender Bildsequenzen geeignet; dies soil hier 
am Beispiel eines Kurzfilms zum graphischen Differenzieren gezeigt werden. Unter Verwendung der 
Hilfsfunktionen LINE, TANGENT, PUNKT und STELLEN defmieren wir die Funktion GD_SCHAR, 
bei der die Parameter / (fiir die zu differenzierende Funktion), a , b (fur die Intervallgrenzen) und n (fur 
die Anzahl der Teilintervalle) gewahlt werden konnen, so dass dem Funktionsbenutzer noch geniigend 
Spielraum bleibt. 

Die hier abgedruckten Bilder konnen DERIVES Leistungen als Filmproduzent naturgemaB nur sehr 
unvollkommen wiedergeben. Der Leser sollte lieber DERIVE starten, die wesentlichen der im Listing 
abgedruckten Zeilen (#3, #5, #7, #9 und #14 sowie #19, #24, #25, #26) eingeben und die Filme star- 
ten. Das Abschreiben ist bekanntlich schnell erledigt, wenn man die F3- und F4-Taste geschickt nutzt. 

Noch ein Tip: Die Ablaufgeschwindigkeit der Filme lasst sich durch den Parameter Accuracy unter 
Plot-Options beeinflussen. 

Und nun viel Vergniigen! 

Das Drehbuch: 


DERIVE is an excellent tool for producing illustrating sequences of graphs. We will demonstrate this 
by an example of a short film illustrating graphic differentiation. By using auxiliary functions LINE, 
TANGENT, PUNKT (point) and STELLEN (places) we define the function GD_SCHAR (gd_family). 
We can choose parameters f ( the given function), a and b (the given interval) and n (the number of 
partial intervals) in order to leave room for investigations. 

Much fun, the script foil wows: 

The other comments are with respect to DERIVE for DOS. I'll come back to them later. Josef 


# 1 : 

# 2 : 

# 3 : 

# 4 : 

#S : 

# 6 : 

# 7 : 

# 8 : 

# 9 : 


Graphisches Differenzieren 

Wi r defimeren die Funktionen LINE, TANGENT, PUNKT, STELLEN und GD_SCHAR : 
LINE(xO, yO, dl, x) := dl-x + yO - dl-xO 


(Hilfsfunktion fiir TANGENT) 


TANGENT(u , x, xO) := LINE 


xO, lira u , lira — u , x 
x-ttO x-»xO dx 


(Term fiir die Tangente in B(xQ;u(xQ)) an den Graphen von u) 

PUNKT(u, xO) := TxO, lim ul 
|_ x-*xO J 

(Punkt auf dem Graphen von u an der Stelle xO) 

(k - l)-(b - a) 


STELLEN (a, b, n) := VECTOR 


■ , k , 1, n + 1 


# 10 : 


(Stellen auf der x-Achse; a: Start, b: Ziel, n: Anzahl der Teilintervalle) 
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#11: 

Bei'spi'el : 



#12: 

STELLENC2, 5, 10) 



#13: 

[2, 2.3, 2.6, 2.9, 3.2, 

3.5, 3.8, 

4.1, 4.4, 

#14: 

GD_SCHARCf, a, b, n) := 

f, VECTOR 

PUNKTCf, 


4.7, 5] 

ELEMENTCSTELLENU, b, n) , k)), TANGENT(f, x, 


f d 11 1 d 1 

ELEMENT(STELLEN(a., b, n) , k)), PUNKT — f, ELEMENT(STELLEN(a., b, n) , k) , k , 1 , n + 1 , — f 

Ldx JJ J dx . 

Zeichnet den Graphen von f ? n aquidistante Beruhrpunkte iiber [a;b] ? die Tangenten an den Graphen von f in diesen Punkten, 
Ableitungspunkte an den Beriihrstellen und zum Schlu(3 den Ableitungsgraphen. 


Plots the graph of f, points on the graph in [a;b], the tangents in this points, the points of the 1 st derivative 
and finally the graph of the first derivative. 

#15 : Ben spiel 1: / Example 1: 

2 

#16: GD_SCHAR(x , -2, 2, 10) 

#17: X fur appnoX engi'bt: 



[-2 

, 4] 

- 4-x - 4 

[-2, 

-4] ' 



[-1.6, 

2.S6] 

- 3.2-x - 2.B6 

[-1.6, 

-3.2] 



[-1.2, 

1.44] 

- 2. 4-x - 1.44 

[-1.2, 

-2.4] 



[-0.8, 

0.64] 

- 1.6-x - 0.64 

[-0.8, 

-1.6] 



[-0.4, 

0.16] 

- O.S-x - 0.16 

[-0.4 , 

-0.8] 


2 







x . 

[0, 

0] 

0 

[0, 

0] 

i 2-x 


[0.4, 

0.16] 

O.S-x - 0.16 

[0.4 , 

0.8] 



[0.8, 

0.64] 

1.6-x - 0.64 

[0.8, 

1.6] 



[1.2, 

1.44] 

2. 4-x - 1.44 

[1.2, 

2.4] 



[1.6, 

2.66] 

3.2-x - 2.56 

[1.6, 

3.2] 



[2, 

4] 

4-x - 4 

[2, 

4] 



Fig. 1 - Fig. 3: Drei Momentaufnahmen des Films GD SCHAR (x A 2, -2, 2, 10). 




Fig. 1: Kurve / graph 


Fig. 2: Kurve mit Tangenten und Ableitungspunkten 

Graph, tangents and points of 1st derivative 
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Fig. 3: Kurve mit Tangenten, Ableitungspunk 
ten und erster Ableitung. 

Graph, tangents, points of 
1st derivative and derivative. 


Drei Schlussszenen three “great movie-finals” 

aus "GD_Schar(x3 , -2.5, 2.5, 50) 11 , 
aus ”GD_Schar(exp(x) , -2.5, 2.5, 20) ", 
aus f, GD_Schar(si n(x) , -x, x, 10)” 






Fig. 6: Vorbereitung einer wichtigen Ableitungsformel / Preparation of an important 1st derivative 


D 


-N-L#6 


Hans-Jiirgen Kayser: Trickfilme mit DERIVE 


p!3 


This was the original contribution from 1992. At this time we had slow computers and we could control 
plotting velocity by changing plotting accuracy. So the plotting process displayed as a movie. Now we 
have very fast computers with fast graphic adapters and we cannot influence plotting velocity. So we 
have to miss the nice movie? 

No, we don't. See the following trick: You must create your own delay function. Plotting shaded areas 
takes some time. So let plot one or more shaded circles far outside of the plot region which is presen- 
ted on the screen. 


2 2 

#24: delay (n) := VECTOR((x - 100) + (y - 100) < RANDOM(l) + k, k, n) 


#25: GD_SCHAR_del(f, a, b, n, t := 0) := 


f, delay tt) , VECTOR 


PUNKTCf, ELEMENTCSTELLENCa, b, n) , k)), 


TANGENTCf, x, ELEMENT(STELLEN(a, b, n) , k)), PUNKT 


I'd 

— f, ELEMENTCSTELLENCa, b, n), k) 
.dx 


dela.y(t) 



d " 

, k, 1, n + 1 

, — f 

J 

dx . 


del ay (n) plots n shaded circles far away. The last parameter of the function controls the plot velocity. 
Omitting this parameter results in a picture - without the movie-effect. t= 1,2 leads to nice animati- 
ons. 

It is worth a try! (Note: Don't simplify the expressions in the Algebra Window, activate Simplify before 
plotting in the 2D-Plt Window.) 


2 


#26: 

GD_SCHAR_del (x , 

-2 


2 


#27: 

GD_SCHAR_del (x , 

-2 


2 


#28 : 

GD_SCHAR_del (x , 

-2 


2 , 10 ) 

2 , 10 , 1 ) 
2 , 10 , 2 ) 


Three stills from the exciting movie GD_SCHAR_del (7x , 0 , 10 , 20 , l) 
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Hyperreelle Zahlen: 

Konstruktion der Zahlenmengen aus den natiirlichen Zahlen: 

a) Aquivalenzklassenbildung 

(welches Kriterium ist jeweils brauchbar?) 

b) Definition der (erweiteren) Rechenarten 
(unabhangig von der Wahl der Reprasentanten moglich?) 

c) Einbettung der bisherigen Zahlen 
(was wird identifiziert?) 

(1, 2, 3, 0, 0, 0, ...) • (0, 0, 0, 4, 4, 4, ...) = (0, 0, 0, 0, 0, 0, ...) 

NULLTEILERKATASTROPHE! ! 

(1,2, 3,0, 0, 0,...) := (0, 0, 0, 0, 0, 0, ...) 

(0, 0, 0, 4, 4, 4,...) :=(4, 4, 4, 4, 4, 4,...) 

( 0 , 1 , 0 , 1 , 0 , 1 , 0 , ...) • ( 1 , 0 , 1 , 0 , 1 , 0 , 1 , ...) = ( 0 , 0 , 0 , 0 , 0 , 0 , 0 , ...) 

NULLTEILERKATASTROPHE! ! 

( 0 , 1 , 0 , 1 , 0 , 1 , 0 ,...) ( 0 , 0 , 0 , 0 , 0 , 0 , 0 , ...) 

(1,0, 1,0, 1,0, 1,...) (1,1, 1,1, 1,1,1,...) 

Grundsatz: 

Es gibt ein Entscheidungsmafi m, das jeder Menge naturlicher Zahlen entweder die Zahl 0 
Oder 1 zuordnet und folgende Eigenschaften hat: 

a) wi(N) — 1 

b) T endlich => «i(T) = 0 

c) T n S = 0 => m(T uS) = wi(T) + m(S) 

Definition: 

Zwei Folgen heiBen Equivalent genau dann, wenn hinreichend viele Fogenglieder 
ubereinstimmen. 

“hinreichend viel”: 

Die Menge der natiirlichen Zahlen, die iibereinstimmende Folgenglieder bezeich- 
nen, hat das EntscheidungsmaB 1. 

Jede Aquivalenzklasse reeller Zahlenfolgen ist eine hyperrelle Zahl. 
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DERIVE lasst sich auch zur (ANALYSIS) Infinitesimalrechnung nutzen. 

Der Standardteil lasst sich wie folgt definieren: 

#1: FFFFFGO := 

ST(fffff, hoi, (ip, 55) := lim lim lim fffff 
#2: 5S->0 (J|J^0 txtx^Q 

Definition der 1. Ableitung FI als Standardanteil des Differentialquotienten 

Definition of 1st derivative FI as standard part of differential quotient 

#3: FGO := 

Different! alquotient 

C lim f) - lim f 
x->a + da x->a 


#4: 

DQ(f, a, 

1 1 
~Tj 


Fl(f, b, 

db) := ST(DQ(f, 

#6: 

3 

FI (x , u 

, du) 

#7: 

2 

3-u 



Integration als Standardanteil ebenfalls for nicht zu komplexe Terme moglich. 

Integration as standard part for not too complex expressions possible. 

Comment: Non-standard analysis was introduced in the early 1960s by the mathematician Abraham 

Robinson. 

Related websites: en.wikipedia.org/wiki/Non-standard_analysis 

en.wikipedia.org/wiki/Non-standard_calculus 
www.math.wisc.edu/~keisler/calc.html 
mathforum.org/dr.math/faq/analysis_hyperreals.html 
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Wunschliste fur Verbesserungen 
Wishlist for Improvements 

Herbert Appel, Schweinfurt 


Wie bereits im DNL#5 berichtet, wurden an der staatlichen Realschule in Forchheim einige Unter- 
richtssequenzen zum Thema „DERIVE im Mathematikunterricht“ durchgefiihrt. Meine Beobachtun- 
gen wahrend dieser Unterrichtsstunden und die abschlieBende Befragung der Schuler brachte einige 
neue Erkenntnisse, die hier aufgezahlt werden sollen. Sie sind als Anregungen bei der Weiterentwick- 
lung von DERIVE und als Forderungen an derartige Programme im Allgemeinen zu verstehen, um 
diese auch fur den Einsatz im Unterricht der Realschule bzw. der Sekundarstufe I geeignet zu machen. 

Obwohl sich die Schuler durchwegs positiv uber das Programm auBerten und sich vorstellen konnten, 
auch in Zukunft mit dem Programm zu arbeiten, ertonten doch auch kritische Stimmen. 

Forderungen beziiglich der Software: 

- Am haufigsten wurden die schwer merkbaren englischen Befehle beanstandet. Weiterhin wurde be- 
klagt, dass das Programm nicht iiber Menus und Maus bedienbar ist. Die Schuler sind von Computem 
(APPLE, AMIGA, ATARI, MS-Windows etc.) anderes (besseres?) gewohnt. 

(The students want to work with menus and to use the mouse.) 

- Schuler tendieren dazu, den Lehrer auf seine Kompetenz hin uberpriifen zu wollen. Nicht benotigte 
Meniipunkte miissen abschaltbar sein! 

(The teacher should be able to switch off menu options which are not needed.) 

- Die Handhabung des Insert/Overwrite-Modus ist bei einigen Meniipunkten inkonsequent. Die zu 
ersetzende Variable wird wird trotz ’’Overwrite” beim Kopieren mit F3 bei Manage Substitute an das 
Zeilenende verschoben. 

(When using the "Overwrite Mode" and copying with F3 to Manage Substitute the expression does not 
overwrite the variable.) 

- Die Handhabung der Befehle Center, Zoom, Scale etc. ist unokonomisch, langwierig und daher un- 
befriedigend (Einsatz der Maus). 

(Working with Center, Zoom, Scale etc. is boring and unsatisfying - could be done better using a mou- 
se.) 


Fiir den Mathematiker stellen viele Eigenheiten des Programms DERIVE keine Schwierigkeiten dar; 
anders ist es bei unseren Schiilern. In theoretischen Uberlegungen wird meist vom idealtypischen, d.h. 
von einem an der Mathematik interessierten, konzentriert arbeitenden, seine Gedanken auf den Pro- 
blemloseprozess ausrichtenden Schuler ausgegangen. An den allgemeinbildenden Schulen treffen wir 
jedoch den „realen“ Schuler an, der haufig bereits bei einfachen Lerninhalten Schwierigkeiten hat. Der 
Teufel steckt wie so oft im Detail. 
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Weitere Forderungen und Anregungen: 

- Die Sichtbarkeit des Multiplikationspunkts sollte einschaltbar sein. (Bsp.: 3-x-y oder 2-2 anstelle 
von 3 x y und 2 2.) 

(It should be possible to switch on the visibility of the multiplication point - 3 x-y or 2-2 instead of 3 x y 
und 2 2) 

- Die Darstellung der Klammem (runde statt der eckigen und geschweiften, bzw. eckige Klammem 
statt der runden) weicht von den allgemeinen Konventionen ab. 

(The use of parentheses and brackets does not follow the usual conventions.) 

- Gleichungen der Form xy/z = wy/v werden nicht automatisch weiter vereinfacht, auch dann nicht, 
wenn y als „positive“ deklariert wurde. 

(Equations of the form xy/z = wy/z will not be simplified automatically.) 

- Es sollte moglich sein, die Variablen (z.B. fur das Verfahren der Parametertransformation 
x ? = 3x + 2 und y ? = x 2 - 5 mit einem " f " versehen zu konnen. 

(It should be possible to write variables with a dash - x',y'.) 

- Die Schreibweise (2x = 4)/2, die beim schrittweisen Umformen der Gleichungen verwendet wird, ist 
vollig unublich. 

(Writing (2x = 4)/2 used when transforming an equation step by step is not common.) 

- Umformungen an trigonometrischen Ausdriicken bleiben fur Schuler durch die Verwendung der 
Signum- und Absolutwertfunktion schwer durchschaubar. 

(Transformations of trigonometric expressions are difficult for pupils because of the use of the SIGN- 
and the ABS-function.) 

- Es ist unbedingt notig vom BogenmaB ins GradmaB umschalten zu konnen, da Schuler kaum mit der 
BogenmaBdarstellung arbeiten. 1} 

(It must be possible to switch from rad to degree, because pupils are not accustomed to work with 
rad .) l) 

- Der in der Geometrie haufig benutzte griechische Buchstabe y kann wegen der bereits vorhandenen 
r-Funktion nicht fur Winkelbezeichnungen verwendet werden. 

(Letter y cannot be used for naming angles, because of the T-function.) 

- Die Angabe der Intervallgrenzen beim Zeichnen eines Graphen sollte handlicher gestaltet werden. 
Insbesondere ware eine Beschriftbarkeit des Graphen wiinschenswert. 

(The declaration of an interval for plotting a graph should be more comfortable, annotations of graphs 
would be great.) 

- Werden mit BUILD zwei Terme voneinander subtrahiert, wobei der zweite ein Bruch mit negativem 
Vorzeichen ist, so folgen in der Darstellung zwei Rechenzeichen aufeinander. 

(Subtracting two expressions using the BUILD-command, we find two following signs, in case of the 
2nd expression being negative.) 


(I tried to translate Mr. Appels wishes for the non-german speaking DUG-members, the editor). 
1} I think you can use the deg-function (e.g. SIN(60 DEG) will display as SIN (60o), the editor) 
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The following article was written by David R. Stoutemyer and was first published by the 

AMERICAN SOCIETY OF MECHANICAL ENGINEERS, NEW YORK, N.Y. 10017 
I am very grateful for the permission to reprint this contribution. 


Mechanical Engineering Applications 
of the DERIVE PC Symbolic Math Program 

David R. Stoutemyer, 

Honolulu, Hawaii 


ABSTRACT 

The Derive® symbolic math program was developed by Soft Ware- 
house, Inc. of Honolulu Hawaii. The program runs on any IBM-PC or NEC 
PC-9801 compatible personal computer having at least 512 kilobytes of 
RAM and MS-DOS version 2.1 or later. Derive features powerful symbolic 
capabilities, 2D and 3D function plots, and numerical methods. These are in- 
tegrated in a menu-driven windowing interface with on-line help, so it is easy 
to learn. The combination of power, low cost, modest hardware require- 
ments, and short learning curve make Derive highly suitable for both engi- 
neering education and research. As illustrations, this article shows, some of 
the ways Derive can be used to assist in the following applications: 

a) random vibration stability, 

b) computing moments of intertia, 

c) robot kinematics, 

d) statics, 

e) automatic control, 

f) computing a radiative heat transfer view factor, 

g) deriving the Hamilton-Jacobi equations of motion for a constrained 
system, and, 

h) deriving the Navier-Stokes equation in paraboloidal coordinates. 

1. INTRODUCTION 

Derive is a menu-driven symbolic math program with on-line help, making it easy to learn 
(References 1, 2, and 7). Derive includes a windowing system that enables you to split the screen into 
windows or overlay them. Figure 1 shows a screen dump in which there is an algebra window contai- 
ning two expressions, a 2D-Plot window showing a parametric plot of the first expression, and a 30- 
Plot window showing a surface plot of the second expression. 

Version 2.0 described here is programmable. Its built-in symbolic capabilities include integrati- 
on, differentiation, solution of algebraic and differential equations, vector algebra and calculus, and 
matrix algebra on matrices having nonnumeric elements. Derive also includes graphic and approxima- 
te numeric capabilities to a wide variety of mechanical engineering problems. Some of these applicati- 
ons use functions defined in auxiliary files distributed with Derive. Others also use functions written 
for this article that are available free on the Derive bulletin board: (217) 337-0926. 
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COMMAND: Algebra Center Eye Focal Grids Hide Length Options rJrrH Quit Uindou 
Zoom 

Enter option 

Center x:G y:0 Length x:10 y:10 Derive 3D-plot 


Figure 1: Derive Screen Dump (Derive 2.0) 



Figure la: Derive Screen Dump (Version 6.10) 
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2. RANDOM VIBRATION STABILITY 

Reference 5 characterizes the stability of a certain system under random vibration according to the 
zeros of the determinant 



-M 

2 

0 

0 

0 

0 

0 

0 

0 

0 


-to 

-P - M 

1 

0 

0 

0 

0 

0 

0 
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0 

- 2-m 

- 2 * p - u 

0 

0 

s 

0 

0 

0 

0 


s 

0 

0 

- 2*b - 

|_J - 2 * w 

0 

0 

0 

0 

0 


0 

0 

0 

1 

-b - M - 

-w 

0 

0 

0 

0 
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0 
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1 
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0 

0 

0 
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- M 
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0 

0 

0 

0 

0 

0 

-to 

0 

-p - M 

1 


0 

0 

0 

0 

0 

0 

t 

-to 

-w 

-p - b - M 

We would like the determinant expressed as 

a polynomial 

in 

//, so we 

declare that ju 

orders ahead 

of any other variables, then simplify this expression. This yields the following product of a fourth 


and a sixth degree polynomial: 


4 3 2 2 2 2 2 

(M + 2 + \\ *(b + p) + n *(b + 3 -bp + p + 2»m + 2-w) + ii-(b + b-(P + 2 -to + 


2 2 

2-w) + 2-p-(ts) + w) - 2 * t) + b -to + b-(p(to + w) - t) + p -w - pt + (to 

2 6 5 4 2 2 

w) )-(m + m * (3 *b + 3*p + 2) + m * (2 *b + 3 -b *(3 -P + 2) + 2 - (p + to + 2-w)) + 

3 2 2 

2 * m *(b * (3 * p +2) + b * (3 *p +3-0) + 2-w) + 2-(p-(to + 3-w) + 2-w)) + 

2 2 2 2 2 
4-m -(b -(P + w) + b-(3-p-(to + w) + 2-w) + 2-w* (P + w)) + 8*M*(b - p - co + 

2 2 
b-w-Cp + m) + 2*p*d)-w) + 16 ■ b * p * ( 0 - w - 4*s ) 

Try deriving this result manually! 
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3. MOMENT OF INERTIA 

Let /u denote a density of some quantity such as mass, which can vary with position. A common 
application of integration is to determine the integral of ju or of its first or second moments along a 
curve or over a plane, curved surface, or volume, ju = 1 is used to compute arclength, area, volume, 
centroid, or the geometric inertia tensor. Applications include center of pressure, center of bouyancy, 
bending of beams, torsion of bars, and rigid body dynamics. File INERTIA.MTH, distributed with 
Derive, contains general purpose functions for computing such quantities in Cartesian, spherical or 
polar cylindrical coordinates. For example, here is one such function definition: 


VOLUME_INERTIA_CARTESIAN(m, X, a, b, y, u, v, z, s, t):: 
b 


• 

t 



* 

'22 



y+z -xy -xz 



2 2 


M 

-xy x+y - y z 



2 2 

, 


-xz -yz x+y 


dz dy dx 


j J u 


Here and throughout this article, displayed expressions are excerpts from files generated by Derive, 
closely matching the appearance of these expressions on the screen or when printed. 

As an example of using this function definition, suppose that we want to compute the volume 
inertia tensor of the intersection of the two homogeneous cylinders: 

2 2 2 

z = r — x 

2 2 2 

z =r - y 

Using ju = 1 for the uniform density, the appropriate expression to simplify is then 

22 22 22 22 

VO L UM E_I N E RTI A_C ARTE SIAN (1 , x, -r, r, y, - /(r -x), J(r -x), z, - /(r - y ), J(_r - y )) 


giving 


112- r 


45 


128 -r 


45 


112- r 


45 


Note that this is a general answer for any r. Such an exact parameterized symbolic solution 
can't be obtained with numerical integration. 


Comment: You cannot find INERTIA.MTH among the utility files distributed with higher versions of 
Derive. It was replaced by ItnegrationApplications.mth. VOLUMEJNERTIA seems to be the same or 
at least very similar: 
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22 22 22 22 
VOLUME_INERTIA(x , -r, r, y, - V(r - x ), V(r - x ), z, - V(r - y ), V(r - y )) 

2 

7 r 

0 0 


2 

8 r 

0 0 

15 

2 

7 r 

0 0 

15 

The results are not the same, but if you multiply the result by the volume of the body, then: 


22 22 22 22 
VOLUME_INERTIA(x , -r, r, y, - V(r - x ), V(r - x ), z, - V(r - y ), V(r - y ))• 

22 22 22 22 
VOLUME(x , -r, r, y, - V(r - x ), V(r - x ), z, - V(r - y ), V(r - y ) , 1) 

5 

112 r 

0 0 


5 

128 r 

0 0 

45 

5 

112 r 

0 0 

45 


4. ROBOT KINEMATICS 

The position and orientation of the tip of a robot arm can be computed by multiplication of matri- 
ces representing joint rotations and translations. For example, reference 3 gives the relevant matrices 
for the Cincinnati T3-756 Robot as 


ml := 


m2 := 


m3 := 


COS (a) 

0 - SIN(of) 

0 ' 


SIN(a) 

0 COS(o<) 

0 


0 

-1 0 


0 


0 

0 0 


1 . 


cos (3) 

- SIN(P) 

0 

a COS(P) 

SIN(P) 

cos(|3) 

0 

a SIN(P) 

0 

0 

1 


0 

0 

0 

0 


1 

COS(0) 

- SIN(0) 

0 

b COS(0) 

SIN(0) 

COS(0) 

0 

b SIN(0) 

0 

0 

1 


0 

0 

0 

0 


1 
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m4 


m5 


m6 


COSCcp) 0 SIN C<p) 0 ' 

SIN Cep) 0 - COS Cep) 0 

0 1 0 0 

0 0 0 1 . 

COS (a) 0 - SIN(a) 0 ' 

SIN(a) 0 COS (a) 0 

0-1 0 0 
0 0 0 1 . 
COS(t) - SIN(t) 0 0 ' 

SIN(t) COS(t) 0 0 

0 0 10 

0 0 0 1 . 


To multiply these matrices we simplify the expression 

ml m2 m3 m4 m5 m6 

giving a matrix having complicated trigonometric expressions as entries. For example, here is the 
3-3 element of this result: 


COS C 3) (SIN CO) SIN((j) COSCcp) + COSO) SIN(cj) SIN(4>)) 

+ SIN ( 3) (COSO) SIN((j) COSO) - SINO) SIN(cj) SIN(4>)) 

To attempt further simplification, we can apply the factor command, giving 

SIN(cr) (COSO) (SIN(0) COSO) + COS(0) SINO)) 

+ SIN ( 3) (COS(0) COSO) - SIN(0) SINO))) 

The parallel structure between the two overall terms in the cofactor of S IN ( a) suggests collec- 
ting products of sinusoids into angle sums in this cofactor. To do so, we can switch trigono-metric 
simplification to collection, highlight the cofactor, then resimplify to obtain the more compact equi- 
valent: 


Manage Tri gononometry Direction: Collect 
SIN(a) SIN(p + (p + 0) 

Comment: Doing the same in Derive 6 we receive another compact equivalent: 

Trigonometry := Collect 
CQS(p +9-0 + 40 CQS((3 + 9 + o + 

2 2 


Again we have a general answer for any angles, rather than a numeric answer for one particular 
set of angles. 
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5. STATICS 

Reference 6 describes a student static problem in which a slender beam of length u and mass m is 
suspended from a point by two cables making angle 0 which each of the two ends of the beam where 
they attach. The problem is to determine the tension f in the remaining cable just after the other cable 
breaks. Let a be the angular acceleration of the beam, and let a be the magnitude of its linear accelera- 
tion. The equations for horizontal, vertical and angular acceleration are 

f COS(0) = m a SIN(0) 


f SIN(0) - m g = m 


a COS(0) 


u (X ^ 


2 


2 

f SIN(0) u m u cx 


2 


12 


Collecting these equations into a vector, then applying the solve command yields the solution 


2 

g-COSO) G.g.SIN(e) g-m-SINO) 

a = , oc = , f = 

2 2 2 
4 - B-COSO) u ■ (4 - B-COSC0) ) 4 - B-COSO) _ 

Note that as with many student problems, the parameters are symbolic rather than numeric. This 
lessens the chance of careless numeric errors masking the instructors’ assessment of the students’ un- 
derstanding, and it makes the answer more general. Computer algebra similarly further lessens the 
chance of careless algebraic errors masking the instructors’ assessment of the students’ modeling ca- 
pabilities. 

Comment: As many of you might know it was David Stoutemyer who was very deep involved in im- 
plementing computer algebra on the TI-92, which has now developed into the Voyage 200. So it is 
nice to see how this problem is solved by the TI-CAS-devices. The solution does not look the same, 
students should find out if the outputs are equivalent or not. 


rat 


. y f:=t y fht t f? y Y ) 

|nigebra|Calc|Qiher|PrgnIQ|Clean Up| j 


1 f ■ cosCe) = a ■ n ■ since) + eql 

f ■ cosCe) = a ■ n ■ since) 

1 f since) - n ■ g = - a ■ cosCe) - + eq2 

- n - (oc ■ u + 2 a - cosCe)) 


f since) - g ■ n = - 


-f- since) -fvua 
2 = ^2~ 


rTTfiY FT? Y fh^ y ft t fs^ t — "i 

ir— |nigebra|Calc|Qiher|PrgnIQ|ciean Up| | 

-f since) 


"N ' U ' OC 
12 


eq3 

-f since) 


■ oc ■ u 


2 12 

■ solueCeql and eq2 and eq3 ? -Cf a oc>) 
_ g ■ cosCe) ^ _ g ■ n since) 


3 (since)) 2 + 1 


and f = - 


3(sinCe)) z + : 





KRYRT 

RAD EXACT 

FUNC H/30 


rat 


Y f:=t t fht t fI Y fg 
nigebr-alCalclQiherlPrgnlQlClean 


i 


-f sinCe)_ - n 
2 


12 


■ + eq3 

~f ■ since) _ ■ a ■ u 

2 12 


■ solyeCeql and eq2 and eq3 ? -Cf a a}) 
g n sinCe) _ , _ 6 ■ g (since)) 2 

:(sinCe)) 2 + 1 [3 (since)) 2 + lj-u 

1 

KRYRT RAD EXACT FUNC H/30 
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6. AUTOMATIC CONTROL 

Reference 4 describes an automatic control example, and here is how it can be done using 
Derive: 

Suppose that we have derived the Laplace domain characteristic expression 

(s + a) (s + b) + k (s + 2 a) (s + 4 b) s 

For stability analysis we can use the expand command to expand and collect similar powers 
of s giving 

3 2 

ks + s (2ak + 4bk + l) + s(a(8bk + l) + b)+ab 

Throughout the analysis we want s to order ahead of the other variables, so we declare that. Next 
we can use the solve command to determine the value of k that makes this expression 0: 

(b + s) ■ (a + s) 
k = 

s ■ (2 ■ a + s) ■ (4 ■ b + s) 

Next to determine the sensitivity of this solution, we can assign the right side to k, then simplify 

k 


d 

s- — k 
ds 

to obtain 

(b + s)'(4«b + s) ■ (a + s)«(2«a + s) 


2 ■ a « (4 ■ b + 2-b-s + s ) + a-s-Cfl-b + b-s + 2-s ) + s -(4-b + 2-b-s + s ) 

Finally, we can use the solve command to determine the values of s that make the sensitivity 0: 

s = - 4-b 

s = -b 


s = -a 



{-4'b -b -2-a ~a> 


MflIH 


RflD EHhlT 


FUNC H/30 
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7. RADIATIVE HEAT TRANSFER 

Radiative heat transfer involves evaluating difficult multiple integrals to compute view factors. 
For example, the view factor from the rectangle 0 < x < a, 0 < y < b, z = 0 of the rectangle x = 0, 
0 <y<b, 0 < z < c is 


c 



0 


Simplifying this expression gives 


2 2 

VCb + c ) -ATAN 

a 

2 2 

^ n/( b + C ) > 

( a ”| fa’l 

b-ATAN ATAN 

L b J L c > 


TT- 1 

2 2 2 
(a - b - c ) - 1 

c 

2 2 

LN(a + b + c 

TT-C TT 

2 2 2 2 2 
: ) (a - b )-LN(a + b ) 

2 2 2 2 
(a - c )-LN(a + c ) 

i 

4 -TT 

2 

-a-c 

2 2 2 

4-TT-a-c 

2 2 2 

4-TT-a-c 

2 

a-LN(a ) (b 

+ c ) -LN(b 

+ c ) b -LN(b ) c-LNCc 

4- 4- 

: ) 


4 -TT-c 4 - 7 T - a - c 4 - tt - a ■ c 4 - tt - a 


Try deriving this result manually, even with the help of integral tables! 

Figure 2 shows a 3D Derive plot of the value of this expression for c = 1, -5 < a < 5, -5 < b < 5. 



Figure 2 
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Figure 2 in Derive 6 
including a slider bar for c 


8. HAMILTON- JACOBI EQUATIONS OF MOTION 

The Hamiltonian approach is one alternative for determining the equations of motion of a con- 
strained mechanical system. Most numerical ordinary-differential-equation solvers require a system of 
first-order equations, and the Hamiltonian approach has the advantage of directly yielding first-order 
rather than second-order differential equations. However, the algebra can be extremely tedious, ma- 
king it a good candidate for computer algebra. Glossing over details: 

a) The potential and kintec energies are expressed in terms of the generalized coordinates and velo- 
cities. 

b) The Lagrangian is computed as the kinetic minus potential energy. 

c) The vector of generalized momenta is defined as the gradient of the Lagrangian with respect to 
the vector of generalized velocities. 

d) These simultaneous algebraic equations are solved for the velocities in terms of the momenta. 

e) The Hamiltonian is then expressed in terms of the generalized coordinates and momenta. 

f) The derivative of the vector of generalized coordinates is then equated to the gradient of the 
Hamiltonian with respect to the vector of momenta. 

g) The derivative of the vector of momenta is then equated to minus the gradient of the Hamiltonian 
with respect to the vector of generalized coordinates. 

File HAMILTON.MTH on the Derive electronic bulletin board automates the process. First to 
accomplish steps a) through e), simplify an expression of the form 

HAMILTONIAN(ke, pe, 0dot, p) 

Here ke is the kinetic energy, pe is the potential energy, 9dot is a vector of variables represen- 
ting the velocities, and p is a vector of variables representing the momenta. 
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9. NAVIER-STOKES EQUATIONS IN PARABOLOIDAL COORDINATES 

Sometimes boundary shapes encourage the use of non-Cartesian coordinates. For example, the 
airflow around a paraboloidal antenna encourages the use of paraboloidal coordinates, making it use- 
ful to derive the governing Navier-Stokes equations in those coordinates. 

The coordinate-free vector form of the Navier-Stokes equation is well known, and Derive has 
built-in functions for determining its differential operators in any orthogonal curvilinear coordinate 
system. Moreover, file COORD.MTH, distributed with Derive, has built-in functions, that determine 
the metrical coefficients hi, h2, and h3 necessary for computing the differential operators in any 
coordinate system. 

Let the transformation from curvilinear coordinates [a, (3, t ] to Cartesian coordinates 

[x, y, z ] be given by 

x := el 
y := e2 
z := e3 


where el, e2 and e3 are expressions depending ona, (3 , and i . Then simplifying an expression of 
the form 

COVARIANT_METRIC_TENSOR( [ el , e2, e3; cx , (3, t]) 
yields the covariant metric tensor g. 

If g is a diagonal matrix, then simplifying an expression of the form 

[a, |3, t ; SQRT_DIAGONAL(g) ] 

yields a geometry matrix whose first row is the curvilinear coordinates and whose second row is the 
corresponding metrical coefficients. 

Such a matrix is an optional second argument for the vector differential operators such as 
LAPLACIAN, GRAD and DIV. If omitted, the default is Cartesian coordinates [x, y, z ] , with unit 
metrical coefficients. 


As an example of this process, to compute the metric tensor for paraboloidal coordinates, we can 
simplify the right side of 





2 2 " 





2 a - [3 


g := CO VAR I ANT _M ETR IC_TE N SO R 

HACOBIAN 


<x-p-T,a-pV(l t), 

i [a. P, t] 


1 


2 

JJ 


giving 


2 2 

a + (3 0 0 

2 2 

a + (3 0 

2 2 
a ■ (3 


2 

1 - T 



0 


0 
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d d d fd ^2 

M VpCa, |3, t, t) — P (a , p , t , t) 4 ■ m ■ — Va (a , p , t , t) 

dp da da l_daj fd f 4 - m 


2 2 


2 

2 2 

2 

\_dv. 

p ' ' 7 

2 2 

3- (ffl + 3 ) 


a-7(a + 

3) 3-(> 

+ P ) 


l. 3 ■ a ■ 

(a + p ) 

Va (a , p , t , t) ^ 

M- j 

"d F 

— V«(a , p 
.dpj 

, T, t) 

f d 


f M 

Vp (a , p, t , t) ^ 


2 2 


2 2 

Ldp 

1 3 J 

2 2 

2 2 


7(a + P ) J 


a + p 



U-(« + 3 ) 

7(a + P ) J 


2 ( d ^2 



d 


2 2 d 


d 

M-(l - t )■ — 

LdT.J 

Va(a, 

3, T, t) 

M-t- — Va(a, p, t, t) 

dT 

M- (a 

+ 2-e ) — VpCa, 

p, t, t) a- 

m — V3 C« . 3, T, t) 
d3 

2 

2 


2 2 


2 2 

2 

2 2 2 

a ■ 

P 


a ■ p 


3 ■ p ■ (« + p ; 

i 

3-0* + 3 ) 

4 2 

jj ■ (4 ■ H + 2 ■ H ■ p 

2 4 

: + p : 

)■¥«(«, p , T , 

2 2 

t) « ■ m ■ (a + 4 ■ p ) ■ 

■Vp(«, P 

, T, t) 



2 

: 2 

2 3 

2 

2 3 




3-a 

■ (a + 

P ) 

3-P-Ca 

+ P ) 




2 

r 

d d 

2 

2 

d 


d 1 

L-P, T - ID-SIGN 


■ M Vt(o< 

, P, T, t) - 3-H-7 (h + 

P ) ■ Vt(< 

S, P, T, t)- Vh( 

;a, p, t, t) - m- 

— Vt((* , 3, T, ti 


k dT da dT dT J 

2 2 2 
3-a -7(« + p )■ |p| 

The other two components are analogous. 

10. CONCLUSIONS 

The above examples range in difficulty from typical student problems through research examples. 
They also range widely throughout mechanical engineering applications. However, eight examples can 
only sparsely sample the space of possible applications. Among the relevant Derive capabilities not 
illustrated here are Fourier series, Taylor series, exact, numeric and series solution of differential equa- 
tions, the use of Bessel and other special functions, approximate numeric integration and root finding. 
Nonetheless, I hope that these examples are varied enough to suggest others in your discipline. I also 
hope that you will send me preprints or reprints of your applications so that I can direct others to your 
work. In fact, I plan to accumulate similar collections of applications in electrical engineering, civil 
engineering, physics and other fields. 
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